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Abstract 

We develop a method for achieving scalable transmission stabilization and switching of N collid¬ 
ing soliton sequences in optical waveguides with broadband delayed Raman response and narrow- 
band nonlinear gain-loss. We show that dynamics of soliton amplitudes in V-sequence transmission 
is described by a generalized V-dimensional predator-prey model. Stability and bifurcation anal¬ 
ysis for the predator-prey model are used to obtain simple conditions on the physical parameters 
for robust transmission stabilization as well as on-off and off-on switching of M out of N soliton 
sequences. Numerical simulations for single-waveguide transmission with a system of N coupled 
nonlinear Schrodinger equations with 2 < V < 4 show excellent agreement with the predator-prey 
model’s predictions and stable propagation over significantly larger distances compared with other 
broadband nonlinear single-waveguide systems. Moreover, stable on-off and off-on switching of mul¬ 
tiple soliton sequences and stable multiple transmission switching events are demonstrated by the 
simulations. We discuss the reasons for the robustness and scalability of transmission stabilization 
and switching in waveguides with broadband delayed Raman response and narrowband nonlinear 
gain-loss, and explain their advantages compared with other broadband nonlinear waveguides. 
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I. INTRODUCTION 


The rates of information transmission throngh broadband optical waveguide links can be 
significantly increased by transmitting many pulse sequences through the same waveguide 
This is achieved by the wavelength-division-multiplexed (WDM) method, where each 
pulse sequence is characterized by the central frequency of its pulses, and is therefore called 


a frequency channel 


^^ications of these WDM or multichannel systems include hber 


waveguides 


0fli, 


optics transmission lines |2|-1^, data transfer between computer processors through silicon 


and multiwavelength lasers 


10 


13| . Since pulses from different frequency 


channels propagate with different group velocities, interchannel pulse collisions are very 
frequent, and can therefore lead to error generation and severe transmission degradation [iHsI, 
On the other hand, the signihcant collision-induced effects can be used for controlling 
the propagation, for tuning of optical pulse parameters, such as amplitude, frequency, and 
phase, and for transmission switching, i.e., the turning on or off of transmission of one 


or more of the pulse sequences 16l42l|. A major advantage of multichannel waveguide 


systems compared with single-channel systems is that the former can simultaneously handle 
a large number of pulses using relatively low pulse energies. One of the most important 
challenges in multichannel transmission concerns the realization of stable scalable control of 
the transmission, which holds for an arbitrary number of frequency channels. In the current 
study we address this challenge, by showing that stable scalable transmission control can 
be achieved in multichannel optical waveguide systems with frequency dependent linear 
gain-loss, broadband delayed Raman response, and narrowband nonlinear gain-loss. 

Interchannel crosstalk, which is the commonly used name for the energy exchange in 
interchannel collisions, is one of the main processes affecting pulse propagation in broadband 
waveguide systems. Two important crosstalk-inducing mechanisms are due to broadband 
delayed Raman response and broadband nonlinear gain-loss. Raman-induced interchannel 
crosstalk is an important impairment in WDM transmission lines employing silica glass 
fibers j^-^, but is also beneficially employed for amplification j^, 1^ . Interchannel 
crosstalk due to cubic loss was shown to be a major factor in error generation in multichannel 
silicon nanowaveguide transmission 33|. Additionally, crosstalk induced by quintic loss can 
lead to transmission degradation and loss of transmission scalability in multichannel optical 
waveguides due to the impact of three-pulse interaction on the crosstalk [l8|, l3^. On the 
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other hand, nonlinear gain-loss crosstalk can be used for achieving energy equalization, 


transmission stabilization, and transmission switching 
In several earlier studies 


17 


16 


20 |. 


2l|, we provided a partial solution to the key problem of 


achieving stable transmission control in multichannel nonlinear waveguide systems, consid¬ 
ering solitons as an example for the optical pulses. Our approach was based on showing 
that the dynamics of soliton amplitudes in iV-sequence transmission can be described by 
Lotka-Volterra (LV) models for N species, where the specihc form of the LV model de¬ 
pends on the nature of the dissipative processes in the waveguide. Stability and bifurcation 
analysis for the steady states of the LV models was used to guide a clever choice of the 
physical parameters, which in turn leads to transmission stabilization, i.e., the amplitudes 


16l-l21|. Furthermore, on- 


of all propagating solitons approach desired predetermined values 
off and off-on transmission switching were demonstrated in two-channel waveguide systems 
with broadband nonlinear gain-loss 19|,l20|. The design of waveguide setups for transmission 
switching was also guided by stability and bifurcation analysis for the steady states of the 
LV models HQ- 


The results of Refs. 


16 


2l| provide the hrst steps toward employing crosstalk induced 


by delayed Raman response or by nonlinear gain-loss for transmission control, stabiliza¬ 
tion, and switching. However, these results are still quite limited, since they do not enable 
scalable transmission stabilization and switching for N pulse sequences with a general N 
value in a single optical waveguide. To explain this, we hrst note that in waveguides with 
broadband delayed Raman response, such as optical hbers, and in waveguides with broad¬ 
band cubic loss, such as silicon waveguides, some or all of the soliton sequences propagate 
in the presence of net linear gain H, 17, 211. This leads to transmission destabilization at 
intermediate distances due to radiative instability and growth of small amplitude waves. As 
a result, the distances along which stable propagation is observed in these single-waveguide 
multichannel systems are relatively small even for small values of the Raman and cubic loss 
coefficients [l7|, 1211]. The radiative instability observed in optical hbers and silicon waveg¬ 
uides can be mitigated by employing waveguides with linear loss, cubic gain, and quintic 


loss, i.e., waveguides with a Ginzburg-Landau (GL) gain-loss prohle jl8l-l2Cll|. However, the 


latter waveguides suher from another serious limitation because of the broadband nature 
of the waveguides nonlinear gain-loss. More specihcally, due to the presence of broadband 
quintic loss, three-pulse interaction gives an important contribution to collision-induced am- 
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plitude shifts 18|, The complex nature of three-pulse interaction in generic three-soliton 


collisions in this case (see Ref. 


34| i leads to a major difficulty in extending the LV model 


for amplitude dynamics from = 2 to a general N value in waveguides with broadband 
nonlinear gain-loss. In the absence of a general iV-dimensional LV model, it is unclear how 
to design setups for stable transmission stabilization and switching in iV-sequence systems 
with N > 2. For this reason, transmission stabilization and switching in waveguides with 


broadband nonlinear gain-loss were so far limited to two-sequence systems 18l-l20[|. 


In view of the limitations of the waveguides studied in Refs. 16l-l2l|. it is important to 


look for new routes for realizing scalable transmission stabilization and switching, which 
work for V-sequence transmission with a general N value. In the current paper we take 
on this task, by studying propagation of N soliton sequences in nonlinear waveguides with 
frequency dependent linear gain-loss, broadband delayed Raman response, and narrowband 
nonlinear gain-loss. Due to the narrowband nature of the nonlinear gain-loss, it affects only 
single-pulse propagation and intrasequence interaction, but does not affect intersequence 
soliton collisions. We show that the combination of Raman-induced amplitude shifts in 
intersequence soliton collisions and single-pulse amplitude shifts due to gain-loss with prop¬ 
erly chosen physical parameter values can be used to realize robust scalable trasmission 
stabilization and switching. For this purpose, we hrst obtain the generalized iV-dimensional 
predator-prey model for amplitude dynamics in an V-sequence system. We then use stabil¬ 
ity and bifurcation analysis for the predator-prey model to obtain simple conditions on the 
values of the physical parameters, which lead to robust transmission stabilization as well as 
on-off and off-on switching of M out of N soliton sequences. The validity of the predator-prey 
model’s predictions is checked by carrying out numerical simulations with the full propaga¬ 
tion model, which consists of a system of N perturbed coupled nonlinear Schrodinger (NLS) 
equations. Our numerical simulations with 2 < V < 4 soliton sequences show excellent 
agreement with the predator-prey model’s predictions and stable propagation over signih- 
cantly larger distances compared with other broadband nonlinear single-waveguide systems. 
Moreover, stable on-off and off-on switching of multiple soliton sequences and stable multiple 
transmission switching events are demonstrated by the simulations. We discuss the reasons 
for the robustness and scalability of transmission stabilization and switching in waveguides 
with broadband delayed Raman response and narrowband nonlinear gain-loss, and explain 
their advantages compared with other broadband nonlinear waveguides. 
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The rest of the paper is organized as follows. In Section m we present the conpled-NLS 
model for propagation of N pulse sequences through waveguides with frequency dependent 
linear gain-loss, broadband delayed Raman response, and narrowband nonlinear gain-loss. 
In addition, we present the corresponding generalized iV-dimensional predator-prey model 
for amplitude dynamics. In Section IIIIl we carry out stability and bifurcation analysis for 
the steady states of the predator-prey model, and use the results to derive conditions on 
the values of the physical parameters for achieving scalable transmission stabilization and 
switching. In Section IIVI we present the results of numerical simulations with the conpled- 
NLS model for transmission stabilization, single switching events, and multiple transmission 
switching. We also analyze these results in comparison with the predictions of the predator- 
prey model. In Section|Vl we discuss the underlying reasons for the robustness and scalability 
of transmission stabilization and switching in waveguides with broadband delayed Raman 
response and narrowband nonlinear gain-loss. Section I VI I is reserved for conclusions. 


II. COUPLED-NLS AND PREDATOR-PREY MODELS 

A. A coupled-NLS model for pulse propagation 

We consider N sequences of optical pulses, each characterized by pulse frequency, propa¬ 
gating in an optical waveguide in the presence of second-order dispersion, Kerr nonlinearity, 
frequency dependent linear gain-loss, broadband delayed Raman response, and narrowband 
nonlinear gain-loss. We assume that the net linear gain-loss is the difference between am- 
pliher gain and waveguide loss and that the frequency differences between all sequences are 
much larger than the spectral width of the pulses. Under these assumptions, the propagation 
is described by the following system of N perturbed coupled-NLS equations: 

N 

+ dlijj + + 6jk) = 

k=l 

N 

igj'ipj/2 + - Sjk) + 'ipkdti'ipj'ipl)] , ( 1 ) 

k=l 

where tjjj is proportional to the envelope of the electric field of the jth sequence, 1 < j < A, 
2 ; is propagation distance, and t is time. In Eq. ([1]), gj is the linear gain-loss coefficient 
for the jth sequence, is the Raman coefficient, and L(|'0jp) is a polynomial of 
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describing the waveguide’s nonlinear gain-loss profile. The values of the gj coefficients are 
determined by the iV-dimensional predator-prey model for amplitude dynamics, by looking 
for steady-state transmission with equal amplitudes for all sequences. The second term on 
the left-hand side of Eq. ([T]) is due to second-order dispersion, while the third and fourth 
terms represent intrasequence and intersequence interaction due to Kerr nonlinearity. The 
Erst term on the right-hand side of Eq. ([1]) is due to linear gain-loss, the second corresponds 
to intrasequence interaction due to nonlinear gain-loss, the third describes Raman-induced 
intrasequence interaction, while the fourth and hfth describe Raman-induced intersequence 
interaction. Since we consider waveguides with broadband delayed Raman response and nar¬ 
rowband nonlinear gain-loss, Raman-induced intersequence interaction is taken into account, 
while intersequence interaction due to nonlinear gain-loss is neglected. The polynomial L 
in Eq. ([1]) can be quite general. In the current paper, we consider two central examples for 
waveguide systems with nonlinear gain-loss: (1) waveguides with a GL gain-loss prohle, (2) 
waveguides with linear gain-loss and cubic loss. The expression for L(|-^jp) for waveguides 
with a GL gain-loss prohle is 

Li{\'ipj?) = , ( 2 ) 

where and eg are the cubic gain and quintic loss coefficients. The expression for 
for waveguides with linear gain-loss and cubic loss is 


L2m") = 


,( 2 ) 


IV’i 


(3) 


where eg is the cubic loss coefficient. We emphasize, however, that our approach can 
be employed to treat a general form of the polynomial L. Note that since some of the 
perturbation terms in the propagation model ([I]) are nonlinear gain or loss terms, the model 
can also be regarded as a coupled system of GL equations. 


The dimensional and c 


laws for NLS solitons jl|. Exactly the same scaling relations were used in our previous 


imensionless physical quantities are related by the standard scaling 


works on soliton propagation in broadband nonlinear waveguide systems 17H21|. In these 


scaling relations, the dimensionless distance 2 ; in Eq. ([T]) is z = X/{2Ld), where X is the 
dimensional distance, = Tq/\/32\ is the dimensional dispersion length, tq is the soliton 
width, and 132 is the second-order dispersion coefficient. The dimensionless retarded time 
is t = t/tq, where r is the retarded time. The solitons spectral width is uq = l/(7r^ro) 
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and the frequency difference between adjacent channels is Az/ = ( 7 rA/ 3 z/o)/ 2 . tfjj = 
where Ej is proportional to the electric field of the jth pulse sequence and Pq is fhe peak 
power. The dimensionless second order dispersion coefficient is d = —1 = / 92 /( 7 Po'ro); where 
7 is the Kerr nonlinearity coefficient. The dimensional linear gain-loss coefficient for the jth 
is related to the dimensionless coefficient via /(yPo). The coefficients 


sequence 


63 ^^ 63 ^^ and 65 are related to the dimensional cubic gain cubic loss pf\ and quintic loss 


Ps, by € 3 ^^ = 2 p 3 ^^/ 7 , e® = 2pf'^/'j, and €5 = 2 P 5 P 0 / 7 , respectively [2^. The dimensionless 
Raman coefficient is = 2tr/tq, where 7 ^ is a dimensional time constant, characterizing 


the waveguide’s delayed Raman response 


The time constant tr can be determined 


from the slope of the Raman gain curve of the waveguide jl. l43|. 


We note that for waveguides with linear gain-loss and cubic loss, some or all of the 
pulse sequences propagate in the presence of net linear gain. This leads to transmission 
destabilization due to radiation emission. The radiative instability can be partially mitigated 
by employing frequency dependent linear gain-loss g{u), z). In this case, the first term on the 
right hand side of Eq. ([T]) is replaced by iE~^{g{u}, z)'ipj)/2, where V’ is the Fourier transform 
of "0 with respect to time, and stands for the inverse Fourier transform. The form of 
g{oJ,z) is chosen such that existence of steady-state transmission with equal amplitudes for 
all sequences is retained, while radiation emission effects are minimized. More specihcally, 
g{u;, z) is equal to a value required to balance amplitude shifts due to nonlinear gain-loss 
and Raman crosstalk, inside a frequency interval of width W centered about the frequency 
of the jth-channel solitons at distance z, /3j{z), and is equal to a negative value gL elsewhere 
35[. Thus, g{oJ,z) is given by: 


f gj if ^j{z) - W/2 < a; < Pj{z) + W/2 for 1 < j < A, 
I gL elsewhere. 


where gi < 0. The width W in Eq. (j4]) satisfies 1 < IK < A/?, where A/3 = /3j+i(0) — /3j(0) 
for 1 < j < A — 1. The values of the gj coefficients are determined by the generalized 
predator-prey model for collision-induced amplitude dynamics, such that amplitude shifts 
due to Raman crosstalk and nonlinear gain-loss are compensated for by the linear gain-loss. 
The values of gL and W are determined by carrying out numerical simulations with Eqs. ©. 
(jS]), and dl]), while looking for the set, which yields the longest stable propagation distance 
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FIG. 1: An example for the frequency dependent linear gain-loss function g{oj,z) of Eq. Q at 
z = 0 in a three-channel system. 


35| . Figured] shows a typical example for the frequency dependent linear gain-loss function 


g{oJ,z) at z = 0 for a three-channel system with gi = 0.0195, g 2 = 0.0267, g^ = 0.0339, 
gL = —0.5, /3i(0) = —15, /52(0) = 0, /53(0) = 15, and W = 10. These parameter values are 
used in the numerical simulations, whose results are shown in Fig. [T] 

The optical pulses in the jth sequence are fundamental solitons of the unperturbed NFS 
equation with central frequency 13j. The envelopes of these solitons are given by z) = 
rjj exp(ixj)sech(xj), where Xj = rjj (t — yj — 2(3jz), Xj = % + ~ Vj) + [v] ~ /^j) ^ind 

r]j, yj, and aj are the soliton amplitude, position, and phase. Due to the large frequency 
differences between the pulse sequences, the solitons undergo a large number of fast inter¬ 
sequence collisions. The energy exchange in the collisions, induced by broadband delayed 
Raman response, can lead to signihcant amplitude shifts and to transmission degradation. 
On the other hand, the combination of Raman-induced amplitude shifts in intersequence 
collisions and single-pulse amplitude shifts due to frequency dependent linear gain-loss and 
narrowband nonlinear gain-loss with properly chosen coefficients can be used to realize ro¬ 
bust scalable transmission stabilization and switching. In the current paper, we demonstrate 
that such stable scalable transmission control can indeed be achieved, even with the simple 
nonlinear gain-loss prohles ([2]) and (|3]). 
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B. A generalized A-dimensional predator-prey model for amplitude dynamics 


The design of waveguide setups for transmission stabilization and switching is based on 
the derivation of LV models for dynamics of soliton amplitudes. For this purpose, we consider 
propagation of N soliton sequences in a waveguide loop, and assume that the frequency 
spacing A(3 between the sequences is a large constant, i.e., A/3 = \l3jj^\{z) — I3j{z)\ S> 1 
for 1 < j < — 1. Similar to Refs. 3, 17|, we can show that amplitude dynamics of 

the N sequences is approximately described by a generalized iV-dimensional predator-prey 
model. The derivation of the predator-prey model is based on the following assumptions. 
(1) The temporal separation T between adjacent solitons in each sequence satisfies: T S> 
1. In addition, the amplitudes are equal for all solitons from the same sequence, but are 
not necessarily equal for solitons from different sequences. This setup corresponds, for 
example, to phase-shift-keyed soliton transmission. (2) As T S> 1, intrasequence interaction 
is exponentially small and is neglected. (3) Delayed Raman response and gain-loss are 
assumed to be weak perturbations. As a result, high-order effects due to radiation emission 
are neglected, in accordance with single-collision analysis. 

Since the pulse sequences are periodic, the amplitudes of all solitons in a given sequence 
undergo the same dynamics. Taking into account collision-induced amplitude shifts due to 
broadband delayed Raman response and single-pulse amplitude changes induced by gain and 
loss, we obtain the following equation for amplitude dynamics of the jth-sequence solitons 
(see Refs. 3, ^ for similar derivations): 


drjj 

dz 


Vj 


N 


9j + F{Vj) + C - j)f{\j - k\)r]k 


k=l 


(5) 


where 1 < j < A^, and C = 4ei?A/3/T. The function F{T]j) on the right hand side of Eq. ([5]) 
is a polynomial in rjj, whose form is determined by the form of L(|'0jp). For Li and L 2 given 
by Eqs. ([2]) and ([3]), we obtain Fi{T]j) = — We^rjj/lb and F 2 {r]j) = /3, 

respectively. The coefficients f{\j — k\) on the right hand side of Eq. ([5]), which describe the 
strength of Raman interaction between jth- and fcth-sequence solitons, are determined by 
the frequency dependence of the Raman gain. For the widely used triangular approximation 
for the Raman gain curve [l, 221, in which the gain is a piecewise linear function of the 


frequency, f{\j — k\) = 1 for 1 < j < N and 1 < k < N 16|. 

In order to demonstrate stable scalable control of soliton propagation, we look for an 
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equilibrium state of the system (|5]) in the form = rj > 0 for 1 < j < iV. Such 
equilibrium state corresponds to steady-state transmission with equal amplitudes for all 
sequences. This requirement leads to: 


N 


g,^-FM-Cr,J2(k-3)f(\3-k\). 


( 6 ) 


Consequently, Eq. (j5]) takes the form 


dVj 

dz 


Vj 


k=l 


N 


nvi) - m+cj2(k- 3)f{\3 - mih - >;) 


k=l 


(7) 


which is a generalized predator-prey model for N species |3a,l37|. Notice that {rj,... ,r]) and 
( 0 ,..., 0 ) are equilibrium states of the model for any positive values of 63 ^^ 65 , 77 , and 

C. 

We point out that the derivation of an iV-dimensional predator-prey model with a gen¬ 
eral N value is enabled by the narrow bandwidth of the waveguide’s nonlinear gain-loss. 
Indeed, due to this property, the gain-loss does not contribute to amplitude shifts in inter¬ 
channel collisions, and therefore, three-pulse interaction can be ignored. This makes the 
extension of the predator-prey model from iV = 2 to a general N value straightforward. As 
a result, extending waveguide setup design from = 2 to a general N value for waveg¬ 
uides with broadband delayed Raman response and narrowband nonlinear gain-loss is also 
straightforward. This situation is very different from the one encountered for waveguides 
with broadband nonlinear gain-loss. In the latter case, interchannel collisions are strongly af¬ 
fected by the nonlinear gain-loss, and three-pulse interaction gives an important contribution 


to the collision-induced amplitude shift 18 


84| . Due to the complex nature of three-pulse 


interaction in gen eric three-soliton collisions in waveguides with broadband nonlinear gain or 


loss (see Ref. 


34|h it is very difficult to extend the LV model for amplitude dynamics from 


= 2 to a generic N value for these waveguides. In the absence of an A^-dimensional LV 
model, it is unclear how to design setups for stable transmission stabilization and switching 
in V-sequence systems with N > 2. As a result, transmission stabilization and switching 


m waveguic 


systems 


es with broadband nonlinear gain-loss have been so far limited to two-sequence 


18|-l20|. 
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III. STABILITY ANALYSIS FOR THE PREDATOR-PREY MODEL ([7]) 


Introduction 


Transmission stabilization and switching are guided by stability analysis of the equilib¬ 
rium states of the predator-prey model ([7|). In particular, in transmission stabilization, we 
require that the equilibrium state {rj,. .. ,7]) is asymptotically stable, so that soliton ampli¬ 
tude values tend to r] with increasing propagation distance for all sequences. Furthermore, 
transmission switching is based on bifurcations of the equilibrium state [r],... ,t]). To ex¬ 
plain this, we denote by rjth the value of the decision level, distinguishing between on and 
off transmission states, and consider transmission switching of M sequences, for example. 
In off-on switching of M sequences, the values of one or more of the physical parameters are 
changed at the switching distance Zg, such that {rj,... ,7]) turns from unstable to asymptoti¬ 
cally stable. As a result, before switching, soliton amplitudes tend to values smaller than 7]th 
in M sequences and to values larger than 7]th in Y — M sequences, while after the switching, 
soliton amplitudes in all N sequences tend to r^, where tj > This means that transmis¬ 
sion of M sequences is turned on at ^ On-off switching of M sequences is realized by 

changing the physical parameters aX z = Zg, such that [r],... ,7]) turns from asymptotically 
stable to unstable, while another equilibrium state with M components smaller than T^th is 
asymptotically stable. Therefore, before switching, soliton amplitudes in all N sequences 
tend to 7], where t] > rjth, while after the switching, soliton amplitudes tend to values smaller 
than rjth in M sequences and to values larger than rjth in Y — M sequences. Thus, transmis¬ 
sion of M sequences is turned off at z > Zg in this case. In both transmission stabilization 
and switching we require that the equilibrium state at the origin is asymptotically stable. 
This requirement is necessary in order to suppress radiative instability due to growth of 


small amplitude waves 18l-l20||. 


The setups of transmission switching that we develop and study in the current paper are 
different from the sing^-pulse switching setups that are commonly considered in nonlinear 
optics (see, e.g.. Ref. [l[ for a description of the latter setups). We therefore point out the 
main differences between the two approaches to switching. First, in the common approach, 
the amplitude value in the off state is close to zero. In contrast, in our approach, the ampli¬ 
tude value in the off state only needs to be smaller than rjth, although the possibility to extend 
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the switching to very small amplitude values does increase switching robustness. Second, 
in the common approach, the switching is based on a single collision or on a small number 
of collisions, and as a result, it often requires high-energy pulses for its implementation. In 
contrast, in our approach, the switching occurs as a result of the cumulative amplitude shift 
in a large number of fast interchannel collisions. Therefore, in this case pulse energies need 
not be high. Third and most important, in the common approach, the switching is carried 
out on a single pulse or on a few pulses. In contrast, in our approach, the switching is carried 
out on all pulses in the waveguide loop (or within a given waveguide span). As a result, the 
switching can be implemented with an arbitrary number of pulses. Because of this property, 
we can refer to transmission switching in our approach as channel switching. Since channel 
switching is carried out for all pulses inside the waveguide loop (or inside a given waveguide 
span), it can be much faster than conventional single-pulse switching. More specihcally, 
channel switching can be faster by a factor oi M x K compared with single-pulse switching, 
where M is the number of channels, whose transmission is switched, and K is the number 
of pulses per channel in the waveguide loop. For example, in a 100-channel system with 
10^ pulses per channel, channel switching can be faster by a factor of 10® compared with 
single-pulse switching. 

Our channel switching approach can be used in any application, in which the same “pro¬ 
cessing” of all pulses within the same channel is required, where here processing can mean 
amplihcation, hltering, routing, signal processing, etc. A simple and widely known exam¬ 
ple for channel switching is provided by transmission recovery, i.e., the amplihcation of a 
sequence of pulses from small amplitudes values below rjth to a desired hnal value above it. 
However, our channel switching approach can actually be used in a much more general and 
sophisticated manner. More specihcally, let pj represent the transmission state of the jth 
channel, i.e., pj = 0 if the jth channel is oh and pj = 1 if the jth channel is on. Then, the 
A^-component vector (pi, ...,Pj, ...,pn), where 1 < j < A^, represents the transmission state 
of the entire A^-channel system. One can then use this A^-component vector to encode infor¬ 
mation about the processing to be carried out on diherent channels in the next “processing 
station” in the transmission line. After this processing has been carried out, the transmission 
state of the system can be switched to a new state, (gi, ...,gAr), which represents the 
type of processing to be carried out in the next processing station. Note that the channel 
switching approach is especially suitable for phase-shift-keyed transmission. Indeed, in this 
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case, the phase is used for encoding the information , an d therefore, no information is lost 
by operating with amplitude values smaller than rjth 


B. Stability analysis for transmission stabilization and off-on switching 


Let us obtain the conditions on the values of the physical parameters for transmission 
stabilization and off-on switching. As explained above, in this case we require that both 
{rj,... ,r]) and the origin are asymptotically stable equilibrium states of the predator-prey 
model (E]). 

We hrst analyze stability of the equilibrium state [rj,... ,r]) in a waveguide with a nar¬ 
rowband GL gain-loss prohle, where F{rij) = Fiirjj). For this purpose, we show that 


N 


( 8 ) 


= ^[rij - 7] + rjln {rj/rij)] , 
i=i 

where t] = {rji,... ,r]j,..., rj^), is a Lyapunov function for Eq. ([7]) 
that > 0 for any r] with r]j > 0 for 1 < j < A^, where equality holds only at the 

equilibrium point. Furthermore, the derivative of Vl along trajectories of Eq. ([7]) satishes: 


Indeed, we observe 


N 


dVi/dz = -(Ibes/lS) + v)iVj - vf x {Vj + - 5k/4) , 

i=i 


(9) 


where k = 65 7 ^ 0. For asymptotic stability, we require dVi/dz < 0. This 

condition is satished in a domain containing {rj,... ,7]) if 0 < k, < Srf /b. Thus, Vl(^) is 
a Lyapunov function for ([7]), and the equilibrium point ( 77 , ..., 77 ) is asymptotically 
stable, if 0 < K < 877^/5 [^. When 0 < k < Arf/b, ( 77 ,... , 77 ) is globally asymptotically 
stable, since in this case, dV^/dz < 0 for any initial condition with nonzero amplitude values. 
When Arf /b < n < 8r]'^/b, dVi/dz < 0 for amplitude values satisfying rjj > {bn/A — 
for 1 < j < A^. Thus, in this case the basin of attraction of {rj, ..., 77 ) can be estimated by 
((5 k/4 — ? 7 ^)^/^, cxo) for 1 < j < A^. For instability, we require dVi/dz > 0 along trajectories 
of ([7]). This inequality is satished in a domain containing ( 77 ,..., 77 ) if k > 8'q‘^/b. Therefore, 
( 77 ,..., 77 ) is unstable for k > 8 ? 7^/5 j^ . 

Consider now the stability properties of the origin for F{T]j) = Fi{rij). Linear stability 
analysis shows that (0 ,..., 0) is asymptotically stable (a stable node) when gj < 0 for 
1 < J < A ^5 be., when all pulse sequences propagate in the presence of net linear loss. To 
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slightly simplify the discussion, we now employ the widely accepted triangular approximation 
'or the Raman gain curve [1, 

n 

16l |. and therefore the net linear gain-loss coefficients take the form 


22|. In this case, f{\j — k\) = 1 ioi 1 < j < N and 1 < A; < iV 


g. = -Fi{ri) - CN{N + 1 ) 77/2 + CNrij. 


( 10 ) 


Since gj is increasing with increasing j, it is sufficient to require < 0. Substituting Eq. 
flTU]) into this inequality, we hnd that the origin is asymptotically stable, provided that 

K> 47/75+ 3C'Ar(iV-l)/(8e577). (11) 

The same simple condition is obtained by showing that VL{r)) = is a Lyapunov 

function for Eq. ([7]). 

Let us discuss the implications of stability analysis for ( 77 ,..., 77 ) and the origin for trans¬ 
mission stabilization and off-on switching. Combining the requirements for asymptotic sta¬ 
bility of both ( 77 ,..., 77 ) and the origin, we expect to observe stable long-distance propagation, 
for which soliton amplitudes in all sequences tend to their steady-state value 77 , provided the 
physical parameters satisfy 

47/75 + 3CN{N - l)/(8e577) <k< 87/75. (12) 


The same condition is required for realizing stable off-on transmission switching. Using 
inequality flT^ . we hnd that the smallest value of 65 , required for transmission stabilization 
and off-on switching, satishes the simple condition 

65 > 15C'A^(iV - 1)/{32t]^). (13) 


As a result, the ratio e^/es should be a small parameter in A^-sequence transmission with 
S> 1. The independence of the stability condition for ( 77 ,... ,77) on N and and the 
simple scaling properties of the stability condition for the origin are essential ingredients in 
enabling robust scalable transmission stabilization and switching. 

Similar stability analysis can be carried out for waveguides with other forms of the non¬ 
linear gain-loss F{r]j) 39|]. Consider the central example of a waveguide with narrowband 
cubic loss, where F{rij) = F 2 {r}j). One can show that in this case VL(r)), given by Eq. (| 8 ]), 
is a Lyapunov function for the predator-prey model ([T]), and that 

N 

dVi/dz = -(4e7/3) + v)ivj - gf < 0, (14) 

1=1 
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for any trajectory with r]j > 0 for 1 < j < iV. Thus, [t], ... ,t]) is globally asymptotically 
stable, regardless of the values of rj, cr, ey, and N. However, linear stability analysis 
shows that the origin is a saddle in this case, i.e., it is unstable. This instability is related 
to the fact that in waveguides with cubic loss, soliton sequences with j values satisfying 
j > {N + l)/2 — 4:ef'^ri/{3CN) propagate under net linear gain, and are thus subject to 
radiative instability. The instability of the origin for uniform waveguides with cubic loss 
makes these waveguides unsuitable for long-distance transmission stabilization. On the other 
hand, the global stability of ( 77 ,..., rj) and its independence on the physical parameters, make 
waveguide spans with narrowband cubic loss very suitable for realizing robust scalable off-on 
switching in hybrid waveguides. To demonstrate this, consider a hybrid waveguide consisting 
of spans with linear gain-loss and cubic loss [F{rij) = ^ 2 ( 77 ^)] and spans with a GL gain-loss 
prohle [F{r]j) = Fi{rij)]. In this case, the global stability of ( 77 ,..., 77 ) for spans with linear 
gain-loss and cubic loss can be used to bring amplitude values close to 77 from small initial 
amplitude values, while the local stability of the origin for spans with a GL gain-loss prohle 
can be employed to stabilize the propagation against radiation emission. 

C. Stability analysis for on-off switching 

We now describe stability analysis for on-off switching in waveguides with a GL gain-loss 
prohle, considering the general case of switching oh of M out of N soliton sequences. As 
explained in Subsection 3.1, in switching oh of M sequences, we require that ( 77 ,... , 77 ) is 
unstable, the origin is asymptotically stable, and another equilibrium state with M com¬ 
ponents smaller than rjih is also asymptotically stable. The requirement for instability of 
[r],... ,ri) and asymptotic stability of the origin leads to the following condition on the 
physical parameter values: 

K > max{ 877 ^/ 5 , 477^/5 -I- 3CN{N — l)/(8e577)}. (15) 

In order to obtain guiding rules for choosing the on-oh transmission switching setups, it 
is useful to consider hrst the case of switching oh of N-1 out of N sequences. Suppose that 
we switch oh the sequences 1 < k < j — 1 and j + 1 < k < N. To realize such switching, we 
require that (0,..., 0, rjgj, 0,..., 0) is a stable equilibrium point of Eq. ([7j). The value of rjsj 
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is determined by the equation 


’f.j - W.Ji - 15ft/(16£5) = 0. 


(16) 


Since the origin is a stable equilibrium point, transmission switching of — 1 sequences 
can be realized by requiring that Eq. (fT6|) has two distinct roots on the positive half of the 
r^j-axis (the largest of which corresponds to r^sj)- This requirement is satished, provided 41| : 


es > 12|9,|/(5 k^). 


(17) 


Assuming that 5^1 < 5^2 < • • • < S'w < 0, we see that the switching off of the — 1 
low-frequency sequences 1 < j < A^ — 1 is the least restrictive, since it can be realized 
with smaller 65 values. For this reason, we choose to adopt the switching setup, in which 
sequences 1 < j < A^ — 1 are switched off. Employing inequality ffT7|l and the triangular- 
approximation-based expression flTOl) for j = N, we hnd that Eq. flT 6 |l has two distinct roots 
on the positive half of the rj^-axis, provided that 

K > {8r]/5) [5k/4 -r]^- 15CN{N - l)/{32e5r])] . (18) 

Therefore, the switching off of sequences 1 < j < A^ — 1 can be realized when conditions 
03 and (fT 8 |) are satished 1421 . 

We now turn to discuss the general case, where transmission of M out of N sequences 
is switched off. Based on the discussion in the previous paragraph, one might expect that 
switching off of M sequences can be most conveniently realized by turning off transmission 
of the low-frequency sequences, 1 < j < M. This expectation is conhrmed by numerical 
solution of the predator-prey model ([7]) and the coupled-NLS model ([T]). For this reason, 
we choose to employ switching off of M sequences, in which transmission in the M lowest 
frequency channels is turned off. Thus, we require that {Q,... ,Q,ris(M+i), ■ ■ ■ ,Vsn) is an 
asymptotically stable equilibrium point of Eq. ([7]). The values of r]s(M+i), ■ ■ ■ ,VsN are 
determined by the following system of equations 

N 

vtj - - l%j/( 16 e 5 ) - 15C'/(16e5)^(A) - j)/(|j - k\)r],k = 0 , (19) 

k=M^l 

where M -|- 1 < j < A^. Employing the triangular approximation for the Raman gain curve 
and using Eq. 03. we can rewrite the system as: 

N 

vtj - - 15 C'/( 16 e 5 ) (fc - j)Vsk -V^ + 

k=M^l 

+ 15CN[{N + l)/2 - j]r//(16e5) = 0. (20) 
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Stability of (0,..., 0, ris(^M+i), ■ ■ ■, Vsn) is determined by calculating the eigenvalues of the 
Jacobian matrix J' at this point. The calculation yields J'jk = 0 for 1 < j < M and j ^ k, 

N 

Jjj = + Ibegr^^/lS - C[N{N + l)r]/2 - ^ krjsk] 

k=M+l 

N 

+C[Nr] - ^ r],k]j for 1 < j < M, (21) 

k=M+l 

Jjk = C{k — j)r]sj for M + 1 < j < N and j ^ k, 

( 22 ) 

and 

N 

= 9j + 4:e^3^Vsj - 166577^3 + CY^{k - j)r]sk for M + 1 < j < (23) 

k=M+l 

Note that the Raman triangular approximation was used to slightly simplify the form of 
Eqs. fl^ - (l23|) . Since J'jk = 0 for 1 < j < M and j 7 ^ k, the hrst M eigenvalues of the 
Jacobian matrix are Xj = Jjj, where the Jjj coefficients are given by Eq. fl2Tl) . Furthermore, 
since Jjj is either monotonically increasing or monotonically decreasing with increasing j, 
to establish stability, it is sufficient to check that either Jmm < 0 or Jn < 0 . To hnd the 
other N — M eigenvalues of the Jacobian matrix, one needs to calculate the determinant of 
the (A^ — M) X (A^ — M) matrix, whose elements are Jjk, where M + 1 < j, A; < A^. The 
latter calculation can also be signihcantly simplihed by noting that for M + 1 < j < A^, 
the diagonal elements are of order 65 , while the off-diagonal elements are of order Nen 
at most. Thus, the leading term in the expression for the determinant is of order . 
The next term in the expansion is the sum oi N — M terms, each of which is of order 
at most. Therefore, the next term in the expansion of the determinant is of 
order {N — at most. Comparing the hrst and second terms, we see that 

the correction term can be neglected, provided that 65 S> We observe that the last 

condition is automatically satished by our on-off transmission switching setup for N ^ 1, 
since stability of the origin requires €5 > N'^eR 3> N^^'^cr [see inequality ([15])] • It follows 
that the other N — M eigenvalues of the Jacobian matrix are well approximated by the 
diagonal elements Jjj for M J- 1 < j < A^. Therefore, for A^ S> 1, stability analysis of 
(0,..., 0, r]s{M+i)) • • •, Vsn) only requires the calculation of A^ — M J- 1 diagonal elements of 
the Jacobian matrix. 
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We point out that the preference for the turning off of transmission of low-frequency 
sequences in on-off switching is a consequence of the nature of the Raman-induced energy 
exchange in soliton collisions. Indeed, Raman crosstalk leads to energy transfer from high- 


frequency solitons to low-frequency ones 
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43l-l48| ■ To compensate for this energy loss 


or gain, high-frequency sequences should be overamplihed while low-frequency sequences 
should be underamplihed compared to mid-frequency sequences 3, 211. As a result, the 
magnitude of the net linear loss is largest for the low-frequency sequences, and therefore, 
on-off switching is easiest to realize for these sequences. It follows that the presence of 
broadband delayed Raman response introduces a preference for turning off the transmission 
of the low-frequency sequences, and by this, enables systematic scalable on-off switching in 
A^-sequence systems. 


IV. NUMERICAL SIMULATIONS WITH THE COUPLED-NLS MODEL 


The predator-prey model ([7]) is based on several simplifying assumptions, which might 
break down with increasing number of channels or at large propagation distances. In partic¬ 
ular, Eq. d?]) neglects the effects of pulse distortion, radiation emission, and intrasequence 
interaction that are incorporated in the full coupled-NLS model ([T]). These effects can lead 
to transmission destabilization and to the breakdown of the predator-prey model descrip¬ 
tion 17H2l|. In addition, during transmission switching, soliton amplitudes can become 


small, and as a result, the magnitude of the linear gain-loss term in Eq. ([T]) might be¬ 
come comparable to the magnitude of the Kerr nonlinearity terms. This can in turn lead 
to the breakdown of the perturbation theory, which is the basis for the derivation of the 
predator-prey model. It is therefore essential to test the validity of the predator-prey model’s 
predictions by carrying out numerical simulations with the full coupled-NLS model ([T]). 

The coupled-NLS system f[T|) is numerically integrated using the split-step method with 
periodic boundary conditions [l|. Due to the usage of periodic boundary conditions, the 
simulations describe pulse propagation in a closed waveguide loop. The initial condition 
for the simulations consists of N periodic sequences of 2K solitons with amplitudes ?7j(0), 
frequencies /5j(0), and zero phases: 

^ ^ Pj(0) exp{i/3j(0)[f - {k + 1/2)T - 5^]} 


k=-K 


cosh{r7j (0) [f — {k + 1/2)T — 5^]} 


(24) 
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where the frequency differences satisfy A/? = /9j+i(0) — /5j(0) 3> 1, for 1 < j < iV — 1. The 
coefficients 6j represent the initial position shift of the jth sequence solitons relative to pulses 
located at {k+l/2)T for —K < k < K—1. To maximize propagation distance in the presence 
of delayed Raman response, we use 6j = {j — 1)T/N for 1 < j < iV. As a concrete example, 
we present the results of numerical simulations for the following set of physical parameters: 
T = 15, A/3 = 15, and A' = 1. In addition, we employ the triangular approximation for the 
Raman gain curve, so that the coefficients f{\j — k\) satisfy f{\j — k\) = 1 for 1 < j,k < N 
la, l2l|]. We emphasize, however, that similar results are obtained with other choices of the 
physical parameter values, satisfying the stability conditions discussed in Section IIIII 

We hrst describe numerical simulations for transmission stabilization in waveguides with 
broadband delayed Raman response and a narrowband GL gain-loss prohle [L(|'0jp) = 
Lid'^jP)] for A = 2, A = 3, and A = 4 sequences. We choose rj = 1 so that the desired 
steady state of the system is (1,..., 1). The Raman coefficient is = 0.0006, while the 
quintic loss coefficient is es = 0.1 for A = 2, 65 = 0.15 for and A = 3, and 65 = 0.25 for 
A = 4. In addition, we choose k = 1.2 and initial amplitudes satisfying ? 7 j( 0 ) > (5k/ 4— 
for 1 < j < A, so that the initial amplitudes belong to the basin of attraction of (1,..., 1). 
The numerical simulations with Eqs. ([T]) and ([ 2 ]) are carried out up to the hnal distances 
= 36110, z /2 = 21320, and Zf^ = 5350, for A = 2, A = 3, and A = 4, respectively. At 
these distances, the onset of transmission destabilization due to radiation emission and pulse 
distortion is observed. The z dependence of soliton amplitudes obtained by the simulations 
is shown in Figs. El^a), El^c), and EJ^e) together with the prediction of the predator-prey 
model (Cl). Figures [ 2 ](b), [ 2 ](d), and [ 2 ](f) show the amplitude dynamics at short distances. 
Figures l3](a), [3](c), and [3](e) show the pulse patterns \'ipj{t,z)\ at a distance z = Zr before 
the onset of transmission instability, where z^ = 36000 for N=2, Zr 2 = 21270 for N=3, 
and Zr 3 = 5300 for N=4. Figures [3](b), [3](d), and [3](f) show the pulse patterns \'ipj(t,z)\ 
at z = Zf, i.e., at the onset of transmission instability. As seen in Fig. d the soliton 
amplitudes tend to the equilibrium value 77 = 1 with increasing distance for A = 2, 3, and 4, 
i.e., the transmission is stable up to the distance z = Zr itl all three cases. The approach to 
the equilibrium state takes place along distances that are much shorter compared with the 
distances along which stable transmission is observed. Furthermore, the agreement between 
the predictions of the predator-prey model and the coupled-NLS simulations is excellent for 
0 < z < Zr- Additionally, as seen in Figs. El^a), dc), andde), the solitons retain their 
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FIG. 2: The z dependence of soliton amplitudes r]j during transmission stabilization in waveguides 
with broadband delayed Raman response and narrowband GL gain-loss for two-sequence [(a) and 
(b)], three-sequence [(c) and (d)], and four-sequence [(e) and (f)] transmission. Graphs (b), (d), 
and (f) show magnified versions of the rjj{z) curves in graphs (a), (c), and (e) at short distances. 
The red circles, green squares, blue up-pointing triangles, and magenta down-pointing triangles 
represent f?i(z), rj 2 {z), ?? 3 (z), and r/ 4 (z), obtained by numerical simulations with Eqs. ([I]) and ([2]). 
The solid brown, dashed gray, dashed-dotted black, and solid-starred orange curves correspond to 
^i('2), V 2 {z), ?/ 3 (^), and r/ 4 (z), obtained by the predator-prey model ([7]). 
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shape at z = Zr despite the large number of intersequence collisions. The distances Zr, along 
which stable propagation is observed, are signihcantly larger compared with those observed 
in other multisequence nonlinear waveguide systems. For example, the value Zr^ = 36000 for 
iV = 2 is larger by a factor of 200 compared with the value obtained in waveguides with linear 
gain and broadband cubic loss Moreover, the stable propagation distances observed in 
the current work for N = 2, N = 3, and iV = 4 are larger by factors of 37.9, 34.3, and 10.6 
compared with the distances obtained in single-waveguide transmission in the presence of 
delayed Raman response and in the absence of nonlinear gain-loss 2l|]. The latter increase 
in the stable transmission distances is quite remarkable, considering the fact that in Ref. 


2 l|, intrasequence frequency-dependent linear gain-loss was employed to further stabilize 


the transmission, whereas in the current work, the gain-loss experienced by each sequence 
is uniform. We also point out that the results of our numerical simulations provide the hrst 
example for stable long-distance propagation of N soliton sequences with > 2 in systems 
described by coupled GL models. 

We note that at the onset of transmission instability, the pulse patterns become distorted, 
where the distortion appears as fast oscillations of I'lpjit, z) \ that are most pronounced at the 
solitons’ tails [see Figs. [3][b), [3](d), and[3][f)]. The degree of pulse distortion is different for 
different pulse sequences. Indeed, for N = 2, the j = I sequence is signihcantly distorted at 
z = Zf^, while no signihcant distortion is observed for the j = 2 sequence. For N = 3, the 
j = 1 sequence is signihcantly distorted, the j = 3 sequence is slightly distorted, while the 
j = 2 sequence is still not distorted at z = Zf^. For = 4, the j = I and j = 4 sequences 
are both signihcantly distorted at z = Zf^, while no signihcant distortion is observed for the 
j = 2 and j = 3 sequences at this distance. 

The distortion of the pulse patterns and the associated transmission destabilization can be 
explained by examination of the Fourier transforms of the pulse patterns z)\. Figure 0] 

shows the Fourier transforms lipjioo, z) \ at z = Zr (before the onset of transmission instability) 
and at z = Zf (at the onset of transmission instability). Figure [5] shows magnihed versions 
of the graphs in Fig. IHfor small \ipj{u,z)\ values. It is seen that the Fourier transforms of 
some of the pulse sequences develop pronounced radiative sidebands at z = Zf. Furthermore, 
the frequencies at which the radiative sidebands attain their maxima are related to the 
central frequencies (3j{z) of the soliton sequences or to the frequency spacing A(3. The 
latter observation indicates that the processes leading to radiative sideband generation are 
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(a) (b) 



FIG. 3: The pulse patterns I'l/jjit, z)\ near the onset of transmission instability for the two-sequence 
[(a) and (b)], three-sequence [(c) and (d)], and four-sequence [(e) and (f)] transmission setups 
considered in Fig. [2l Graphs (a), (c), and (e) show |'0j(t,z)| before the onset of instability, while 
graphs (b), (d), and (f) show |'0j(t,z)| at the onset of instability. The solid red curve, dashed- 
dotted green curve, blue crosses, and dashed magenta curve represent |V^j(t, z)| with j = 1,2,3,4, 
obtained by numerical solution of Eqs. © and (|2]) . The propagation distances are z = = 36000 

(a), z = z/j = 36110 (b), z = z^j = 21270 (c), z = z/j = 21320 (d), z = z^g = 5300 (e), and 
z = zf^= 5350 (f). 
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resonant in nature (see also Refs. HQ)- 

Consider first the Fourier transforms of the pulse patterns for N = 2. As seen in Figs. 
111(b) and[H(b), in this case the j = I sequence develops radiative sidebands at frequencies 
= 17.18 and = 34.76 at z = Zf^. In contrast, no signihcant sidebands are 
observed for the j = 2 sequence at this distance. These findings explain the significant 
pulse pattern distortion of the j = I sequence and the absence of pulse pattern distortion 
for the j = 2 sequence at z = Zf^. In addition, the radiative sideband frequencies satisfy 
the simple relations: — /d 2 ( 2 ^ri) ~ 29.3 ~ 2A/3 and ~ For N = 3, the 

j = I sequence develops signihcant sidebands at frequencies = 0.0 and = 44.4, 
the j = 3 sequence develops a weak sideband at frequency oos = —31.42, and the j = 2 
sequence does not have any signihcant sidebands at ^ = Zf^ [see Figs. 111(d) and|H(d)]. These 
results coincide with the signihcant pulse pattern distortion of the j = 1 sequence, the weak 
pulse pattern distortion of the j = 3 sequence, and the absence of pulse pattern distortion 
for the j = 2 sequence at z = zj^. Additionally, the sideband frequencies satisfy the simple 
relations: ~ P 3 {zr 2 ), ~ 3A/3, and ~ f3i{zr2)- For N = 4, the j = 1 and 

j = 4 sequences develop signihcant sidebands, while no signihcant sidebands are observed 
for the j = 2 and j = 3 sequences at z = Zf^ [see Figs. 111(f) and|5l(f)]. These hndings explain 
the signihcant pulse pattern distortion of the j = I and j = 4 sequences and the absence 
of signihcant pulse pattern distortion for the j = 2 and j = 3 sequences at z = Zf^. The 
sideband frequencies of the j = 1 sequence satisfy the relations: = 17.17 Psizra), and 

= 34.77 ~ 2uji^^\ Note that the values of and for A = 4 are very close to 
the values found for N = 2. Finally, the sideband frequencies of the j = 4 sequence satisfy 
the relations: = —27.65 ~ (3i{zr3), and = 34.35 ~ 2uji^^K 

We now turn to describe numerical simulations for a single transmission switching event 
in waveguides with broadband delayed Raman response and a narrowband GL gain-loss 
prohle. As described in Section m on-oh switching of M out of N pulse sequences at a 
distance z = Zg is realized by changing the value of one or more of the physical parameters, 
such that the steady state {rj,... ,t]) turns from asymptotically stable to unstable, while 
another steady state at (0,..., 0, 7]s(^m+i), ■ ■ ■, Vsn) is asymptotically stable. We denote the 
on-oh switching setups by A1-A2, where Al and A2 denote the sets of physical parameters 
used at 0 < z < Zg and z > Zg, respectively. 

Oh-on switching of M out of N soliton sequences at z = Zg is realized by changing the 
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(a) (b) 



FIG. 4: The Fourier transforms of the pulse patterns \iljj[uj,z)\ near the onset of transmission 
instability for the two-sequence [(a) and (b)], three-sequence [(c) and (d)], and four-sequence [(e) 
and (f)] transmission setups considered in Figs. [5]and[31 Graphs (a), (c), and (e) show z)\ 

before the onset of instability, while graphs (b), (d), and (f) show z)\ at the onset of in¬ 

stability. The red circles, green squares, blue up-pointing triangles, and magenta down-pointing 
triangles represent |'0j(w, z)\ with j = 1, 2,3,4, obtained by numerical solution of Eqs. ([H) and ([2]). 
The propagation distances are z = Zr^ = 36000 (a), z = Zf^ = 36110 (b), z = Zr 2 = 21270 (c), 
z = Zf^ = 21320 (d), z = Zr^ = 5300 (e), and z = Zf^ = 5350 (f). 
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(d) 




FIG. 5: Magnified versions of the graphs in Fig. S]for small z)\ values. The symbols and 

distances are the same as in Fig. SI 


physical parameter values such that {rj,... ,7]) turns from unstable to asymptotically stable. 
As explained in Section IIIIl to achieve stable long-distance transmission after the switching, 
one needs to require that the origin is an asymptotically stable steady state as well. Under 
this requirement, k must satisfy inequality flT^ . and as a result, the basin of attraction of 
(?],..., 77 ) is limited to {{5n/A — 77 ^)^/^, 00 ) for 1 < j < N. This leads to limitations on 
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the turning on of the M sequences, especially for M > 2 and N > 3. To overcome this 
difficulty, we consider a hybrid waveguide consisting of a span with a GL gain-loss prohle, a 
span with linear gain-loss and cubic loss, and a second span with a GL gain-loss prohle. The 
introduction of the intermediate waveguide span with linear gain-loss and cubic loss enables 
the turning on of the M sequences from low amplitude values due to the global stability of 
the steady state {rj,... ,7]) for the corresponding predator-prey model. However, due to the 
presence of linear gain and the instability of the origin for the same predator-prey model, 
propagation in the waveguide span with linear gain-loss and cubic loss is unstable against 
emission of small amplitude waves. For this reason, we introduce the frequency dependent 
linear gain-loss g{u,z) of Eq. (jTj) when simulating propagation in the second span. More 
importantly, propagation in the second span with a GL gain-loss prohle leads to mitigation of 
radiative instability due to the presence of linear loss in all channels for this waveguide span. 
This enables stable long-distance propagation of the N soliton sequences after the switching. 
We denote the oh-on switching setups by A2-B-A1, where A2, B, and Al denote the sets of 
physical parameters used in the hrst, second, and third spans of the hybrid waveguide. The 
hrst span is located at [0,Zsi), the second at [^51,^52), and the third at [zs 2 ,Zf], where Zf is 
the hnal propagation distance. Thus, oh-on switching of the M soliton sequences occurs at 
z > Zsi, while hnal transmission stabilization takes place at ^ > Zs 2 - 

We present here the results of numerical simulations for on-oh and oh-on switching of 
two and three soliton sequences in four-sequence transmission. As discussed in the preceding 
paragraphs, on-oh switching setups are denoted by A1-A2 and oh-on switching setups are 
denoted by A2-B-A1. The following values of the physical parameters are used. The Raman 
coefficient is = 0.0006, which is the same value used in transmission stabilization. The 
other parameter values used in setup Al in both on-oh and oh-on switching are 65 = 0.1, 
K = 1.2, and r] = 1. The parameter values used in setup A2 in on-oh switching are €5 = 0.04, 
K = 1.8, and rj = 1.05 for M = 2 and 65 = 0.04, k = 2, and 77 = 1.1 for M = 3. The on-oh 
switching distance is Zg = 250 for both M = 2 and M = 3. The parameter values used in 
setup A2 in oh-on switching are €5 = 0.032, k = 2.2, and rj = 1.1 for M = 2 and 65 = 0.02, 
K = 2.8, and r] = 1.3 for M = 3. The parameter values used in oh-on switching in setup 
B are eg = 0.02 and rj = 1 for both M = 2 and M = 3. To suppress radiative instability 
during propagation in waveguide spans with linear gain-loss and cubic loss (setup B), the 
frequency dependent linear gain-loss g{u;, z) of Eq. (jl]) with W = 10 and g^ = —0.5 is 
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employed. The switching and hnal stabilization distances in off-on transmission switching 
are Zgi = 30 and Zs 2 = 80 for M = 2, and Zgi = 30 and Zs 2 = 90 for M = 3. We point 
out that similar results were obtained with other choices of the physical parameter values, 
satisfying the stability conditions discussed in Section IIIII 

The results of numerical simulations with Eqs. ([T]) and ([2]) for on-off switching of two and 
three soliton sequences in four-sequence transmission in setup A 1 -A 2 are shown in Figs. | 6 ]^a) 
and[ 6 ](b). The results of simulations with Eqs. ([I])-(I1]) for off-on switching of two and three 
sequences in four-sequence transmission in setup A2-B-A1 are shown in Figs. [ 6 ]^c) and[ 6 ](d). 
A comparison with the predictions of the predator-prey model ([7]) is also presented. The 
agreement between the coupled-NLS simulations and the LV model’s predictions is excellent 
in all four cases. More specihcally, in on-off transmission of M sequences with M = 2 and 
M = 3, the amplitudes of the solitons in the M lowest frequency channels tend to zero, 
while the amplitudes of the solitons in the N — M high frequency channels tend to new 
values T]sj, where M -|- 1 < j < A^. The values of the new amplitudes are rjss = 1.2499 
and t]sa = 1.2878 in on-off switching of two sequences, and 77^4 = 1.3640 in on-off switching 
of three sequences. As can be seen from Figs. ma) and[ 6 ](b), these values are in excellent 
agreement with the predictions of the predator-prey model ([7]). In off-on switching of M 
soliton sequences, the amplitudes of the solitons in the M low frequency channels tend to 
zero for z < Zgi, while the amplitudes of the solitons in the N — M high frequency channels 
increase with z for z < Zgi- After the switching, i.e., for distances z > Zgi, the amplitudes of 
the solitons in the M low frequency channels increase to the steady-state value of 1, while 
the amplitudes of the solitons in the N — M high frequency channels decrease and tend 
to 1, in full agreement with the predator-prey model’s predictions. Note that very good 
agreement between the coupled-NLS and predator-prey models is observed even when some 
of the soliton amplitudes are small, i.e., even outside of the perturbative regime, where the 
predator-prey model is expected to hold. The results of the simulations presented in Fig. 
m and similar results obtained with other sets of the physical parameters demonstrate that 
it is indeed possible to realize stable scalable on-off and off-on transmission switching in 
the waveguide setups considered in the current study. Furthermore, the simulations conhrm 
that design of the switching setups can be guided by stability and bifurcation analysis for 
the steady states of the predator-prey model ([7]). We point out that the off-on switching 
setups can also be employed in broadband transmission recovery, that is, in the simultaneous 
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FIG. 6: The z dependence of soliton amplitudes r]j during single transmission switching events in 
four-sequence transmission in waveguides with broadband delayed Raman response and narrow- 
band nonlinear gain-loss, (a) and (b) show on-off switching of two and three soliton sequences in 
waveguides with a GL gain-loss profile in setup A1-A2, while (c) and (d) show off-on switching of 
two and three sequences in hybrid waveguides in setup A2-B-A1. The red circles, green squares, 
blue up-pointing triangles, and magenta down-pointing triangles represent 'r] 2 {z), rjsiz), and 

ri 4 {z), obtained by numerical simulations with Eqs. ([I|) and ([2]) in (a) and (b), and with Eqs. 
([I])-([l]) in (c) and (d). The solid brown, dashed gray, dashed-dotted black, and solid-starred orange 
curves correspond to t]i{z), r]2{z), fl^{z), and rj4^{z), obtained by the predator-prey model (j7jl. 

amplification of multiple soliton sequences, which experienced significant energy decay, to a 
desired steady-state amplitude value. 

As discussed in Section imi an important application of the switching setups considered 
in our paper is for realizing efficient signal processing in multichannel transmission. In such 
application, the amplitude values rjj are used to encode information about the type of signal 
processing to be carried out in the next processing station. As a result, the pulse sequences 
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typically undergo multiple switching events and it is important to show that this can be re¬ 
alized in a stable manner. We therefore turn to discuss the results of numerical simulations 
with the coupled-NLS model ([T]) for multiple switching events. As a specihc example, we 
consider multiple switching in a three-channel system in the hybrid waveguide setup Al- 
(A2-B-Al)-...-(A2-B-Al), where (A2-B-A1) repeats six times. Thus, in this case the soliton 
sequences hrst experience transmission stabilization in waveguide setup Al, and then un¬ 
dergo six successive off-on switching events in waveguide setup A2-B-A1. The parameter 
values are chosen such that during on switching, amplitude values in all sequences tend to 1 , 
while during off switching, transmission of a single seqnence (the lowest-frequency sequence 
j = 1) is tnrned off. We emphasize, however, that similar resnlts are obtained with other 
nnmbers of channels and in other switching scenarios. For the example presented here, dur¬ 
ing on switching stabilization in waveguide setup Al, eR = 0.0006, 65 = 0.15, n = 1 . 2 , and 
T] = 1 are used, while during off switching in waveguide setup A 2 , sr = 0.0024, 65 = 0.024, 
/t = 2 . 2 , and rj = 1.15 are used. Note that the higher value of er in waveguide setup A 2 
is required for realizing a faster on-off transmission switching, that is, for decreasing the 

distance along which the off switching takes place. Additionally, dnring on switching in 

f2) 

wavegnide setnp B, cr = 0.0006, 63 = 0.02, and rj = 1 are chosen. To suppress radiative 
instability dnring propagation in waveguide spans with linear gain-loss and cubic loss (setup 
B), the frequency dependent linear gain-loss g{u, z) of Eq. (jl]) with W = 10 and qr = —0.5 
is employed. In the simulations, transmission of the j = 1 sequence is turned off in setnp 
A2 at distances z^m+i = 700(m -|- 1) for 0 < m < 5. Transmission of seqnence j = 1 is 
turned on in waveguide setup B at z^rn +2 = 50-|-700(m-|-l) for 0 < m < 5, and transmission 
stabilization in setnp Al starts at Zq = 0 and at z^m = 100 -|- 700m for 1 < m < 6. The 
hnal propagation distance is Zf = 5000. Thus, the waveguide spans are [0,700), [700,750), 
[750,800), [800,1400), ..., [4200,4250), [4250,4300), and [4300,5000]. 

The resnlts of nnmerical simnlations with Eqs. 0-0 for mnltiple transmission switching 
events are shown in Fig. [7] along with the predictions of the predator-prey model ((71). The 
agreement between the conpled-NLS simulations and the predator-prey model’s predictions 
is excellent thronghout the propagation. Fnrthermore, as seen in Fig. [71(b), no pulse dis¬ 
tortion is observed at the hnal propagation distance Zf. Note that the minimal valnes of 
r]i{z), which are attained prior to the start of the on switch, are r]i{z 3 rn+ 2 ) = 0.33. As a 
result, the value of the decision level rjth for distinguishing between on and off transmission 
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states can be set as low as rjth = 0.35, which is signihcantly lower than the valne rjth = 0.65 


obtained in Ref. 


19[ | for transmission in a two-channel wavegnide system with a broadband 


GL gain-loss prohle. The resnlts presented in Fig. [7] along with resnlts of nnmerical sim- 
nlations with other sets of the physical parameter valnes demonstrate that stable multiple 
transmission switching events can indeed be realized over a wide range of amplitude val¬ 
ues, using waveguides with broadband delayed Raman response and narrowband nonlinear 
gain-loss. 


V. DISCUSSION 


Let us discuss the reasons for the robustness and scalability of transmission stabilization 
and switchng in waveguides with broadband delayed Raman response and narrowband non¬ 
linear gain-loss. The scalability and robustness of transmission control can be attributed 
to the following properties of these waveguides. (1) The asymptotic stability of the steady 
state {rj,... , 1 ]) for waveguides with a GL gain-loss prohle, which is independent of N and 
en, is key to realizing scalable transmission stabilization and switching. (2) The presence 
of net linear loss in all frequency channels for waveguides with a GL gain-loss prohle leads 
to mitigation of radiative instability. (3) Due to the narrow bandwidth of the nonlinear 
gain-loss, three-pulse interaction does not contribute to collision-induced amplitude shifts. 
As a result, the extension of the predator-prey model from = 2 to a general N value 
is straightforward. This also makes the extension of waveguide setup design from N = 2 
to a general N value straightforward. In contrast, in waveguides with broadband nonlin¬ 
ear gain-loss, thr ee-p ulse interaction gives an important contribution to collision-induced 


amplitude shifts 


18 


34| . Due to the complex nature of three-pulse interaction in generic 


three-soliton collisions in waveguides with broadband nonlinear gain or loss (see Ref. 34|L 
it is very difficult to extend the LV model for amplitude dynamics from = 2 to a generic N 
value in these waveguides. In the absence of an A^-dimensional LV model, it is unclear how 
to design setups for stable transmission stabilization and switching in V-sequence waveg¬ 
uide systems with broadband nonlinear gain-loss. As a result, transmission stabilization 
and switching in waveguides with broadband nonlinear gain-loss have been so far limited 


to two-sequence systems 


18 


20| . (4) The Raman-induced energy transfer in soliton colli¬ 


sions from high-frequency solitons to low-frequency solitons is an important ingerdient in 
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FIG. 7: (a) The z dependence of soliton amplitudes rjj in multiple transmission switching with 
three sequences {N = 3) in hybrid waveguide setup Al-(A2-B-Al)-...-(A2-B-Al), where (A2-B-A1) 
repeats six times. In this case, the soliton sequences undergo transmission stabilization followed 
by six successive off-on switching events. The red circles, green squares, and blue up-pointing 
triangles represent r]i{z), rj 2 {z), and r] 3 {z), obtained by numerical simulations with Eqs. ({T ]) -(| ^ . 
The solid brown, dashed gray, and dashed-dotted black curves correspond to r]i{z), rj 2 {z), and 
7 / 3 ( 2 ;), obtained by the predator-prey model ([7]). (b) The pulse patterns at the final 

distance z = Zf = 5000, as obtained by numerical solution of Eqs. ([I])-(I1]). The solid red curve, 
dashed-dotted green curve, and blue crosses represent \'4)j{t,z)\ with j = 1,2,3. 
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the realization of scalable on-off switching. Indeed, to compensate for the Raman-induced 
energy loss or gain in the collisions, high-frequency sequences should be overamplihed while 
low-frequency sequences should be underamplihed compared to mid-frequency sequences. 
As a result, the magnitude of the net linear loss is largest for the low-frequency sequences, 
and therefore, on-off switching is easiest to realize for these sequences. Thus, the presence of 
broadband delayed Raman response introduces a preference for turning off the transmission 
of the low-frequency sequences, and by this enables systematic scalable on-off switching. 
(5) The global asymptotic stability of the steady state {rj,... ,7]) for waveguide spans with 
linear gain-loss and cubic loss, which is independent of all physical parameters, is important 
for realizing robust scalable off-on switching in hybrid waveguides. These five waveguide 
properties are explained by stability and bifurcation analysis for the steady states of the 
generalized A^-dimensional predator-prey model for amplitude dynamics. Thus, the analysis 
of the predator-prey model is essential to the design of waveguide setups leading to stable 
scalable control of soliton-based multichannel transmission. 

Note that waveguide setups with narrowband cubic gain and quinitc loss or with nar¬ 
rowband cubic loss can be realized by employing fast saturable absorbers with a bandwidth 
Auci satisfying: z/q ^ An. That is, in these waveguide systems, the bandwidth of 

the saturable absorber is larger than the spectral width of the optical pulses but smaller than 
the frequency spacing between adjacent frequency channels. Waveguide systems containing 
a fast saturable absorber with a hnite spectral width, which is much larger than the spectral 
width of the optical pulses, have been studied extensively in the context of mode-locked 


lasers; see, for example. Refs. 50|-|52| and references therein 


VI. CONCLUSIONS 

We developed a method for achieving stable scalable control of propagation of multi¬ 
ple soliton sequences in broadband optical waveguide systems. The method is based on 
employing nonlinear waveguides with broadband delayed Raman response, linear gain-loss, 
and narrowband nonlinear gain-loss. We showed that the combination of Raman-induced 
amplitude shifts in interchannel collisions and single-pulse amplitude shifts due to linear and 
nonlinear gain-loss with properly chosen physical parameter values can be used to realize 
robust scalable transmission stabilization and switching. For this purpose, we first showed 
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that the dynamics of soliton amplitudes in an iV-sequence transmission system can be de¬ 
scribed by a generalized iV-dimensional predator-prey model. We then carried out stability 
and bifurcation analysis for the steady states of the predator-prey model for two central 
cases of the gain-loss: (1) a GL gain-loss prohle, (2) linear gain-loss and cubic loss. The 
stability and bifurcation analysis was then used to develop waveguide setups that lead to 
robust transmission stabilization as well as on-off and off-on switching of M out of N soliton 
sequences. 

For waveguides with a GL gain-loss prohle, we obtained the Lyapunov function Vl(^) for 
the predator-prey model and used it to derive simple conditions for asymptotic stability and 
instability of the steady state with equal amplitudes for all sequences {rj,... ,r]). These con¬ 
ditions are independent of the number of channels N and the value of the Raman coefficient 
epi, which is essential for the realization of scalable transmission stabilization and switch¬ 
ing. We also found that the steady state at the origin is asymptotically stable, provided 
all the linear gain-loss coefficients are negative. Gombining the requirements for asymptotic 
stability of both [rj,... , 7 ]) and the origin, we showed that the smallest value of the quin- 
tic loss coefficient es required for robust transmission stabilization and off-on switching for 

S> 1 scales as 65 ~ N‘^€-r- The realization of on-off switching requires stability analysis of 
steady states, for which M components are equal to zero. We hrst gave a simple argument, 
showing that switching off of M sequences is most conveniently realized by turning off the 
transmission of the low-frequency sequences 1 < j < M. We therefore focused attention 
on the steady state ( 0 ,..., 0 , ris(M+i), ■ ■ ■, Vsn) and showed that for iV S> 1, stability of this 
steady state can be established by calculating only N — M + 1 diagonal elements of the 
corresponding Jacobian matrix. 

Stability analysis for the predator-prey model, describing amplitude dynamics in waveg¬ 
uides with linear gain-loss and cubic loss, was carried out in a similar manner. More specif¬ 
ically, we found that the same Vl(j)) that was used for the predator-prey model with a GL 
gain-loss prohle is a Lyapunov function for the predator-prey model with linear gain-loss 
and cubic loss |39[. Moreover, we used this Lyapunov function to show that {rj,... , 1 ]) is 
globally asymptotically stable, regardless of the values of all physical parameters. However, 
linear stability analysis showed that the origin is unstable in this case. The latter instability 
eventually leads to growth of small amplitude waves, and thus makes waveguides with lin¬ 
ear gain-loss and cubic loss unsuitable for long-distance transmission stabilization. On the 
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other hand, the global asymptotic stability of ( 77 ,, rj) means that waveguide spans with 
linear gain-loss and cubic loss can be used in hybrid waveguides for realizing robust off-on 
transmission switching. 

The predictions of the generalized predator-prey model for scalable transmission stabi¬ 
lization and switching were tested by numerical simulations with a perturbed coupled-NLS 
model, which takes into account broadband delayed Raman response and a narrowband GL 
gain-loss prohle. The coupled-NLS simulations for transmission stabilization were carried 
out with 2 < iV < 4 soliton sequences. The simulations showed stable propagation and 
excellent agreement with the predictions of the predator-prey model over signihcantly larger 
distances compared with those obtained in earlier works with other waveguide setups. More 
specihcally, the stable propagation distances obtained for two-, three-, and four-sequence 
transmission were larger by factors of 37.9, 34.3, and 10.6, respectively, compared with 
the distances obtained in single-waveguide transmission in the presence of delayed Raman 
response and in the absence of nonlinear gain-loss 2l|]. Furthermore, the distance along 
which stable transmission was observed in a two-channel system was larger by a factor of 
200 compared with the distance achieved in waveguides with linear gain and broadband 


cubic loss 


17l |. The enhanced stability of iV-channel transmission through waveguides with 


broadband delayed Raman response and narrowband GL gain-loss prohle was explained in 
the Discussion. 

We demonstrated single and multiple transmission switching events of M out of N pulse 
sequences by carrying out numerical simulations with the coupled-NLS model that was de¬ 
scribed in the preceding paragraph. As examples, we presented the results of the simulations 
for the following setups: (a) single on-off and off-on switching events of two and three soli- 
ton sequences in four-sequence transmission; (b) six switching events in a three-sequence 
system, in which transmission of one soliton sequence was switched off and then on six con¬ 
secutive times. The results of the coupled-NLS simulations were in excellent agreement with 
the predictions of the predator-prey model for both single and multiple switching events. 
Furthermore, the agreement was observed even when amplitude values were small for some 
soliton sequences, i.e., even outside of the regime where the predator-prey model’s descrip¬ 
tion was expected to hold. Based on these results and results of simulations with other sets of 
the physical parameter values, we concluded that stable scalable transmission switching can 
indeed be realized in waveguides with broadband delayed Raman response and narrowband 
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nonlinear gain-loss. 
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